AD-A227  851 


■'»  r 


^OPV  Si 


Julius  Silver  Institute  of  Biomedical  Engineering 


Technion,  Israel  Institute  ofTechnology,  Haifa,  Israel  32000 


GRANT  AF0SR-89-0054 


CEREBRAL-BODY  PERFUSION  MODEL 


S.  Sorek1,  J.  Bear2,  and  ML  Feinsod3 
in  Collaboration  with 

K.  Allen4,  L.  Bunt5  and  S.  Ben-Haim6 


July  1990 

SCIENTIFIC  REPORT  NO.  4 


Approved  for  public  release,  distribution  unlimited. 


Prepared  for 

The  United  States  Air  Force,  Air  Force  Office  of  Scientific  Research 

and 


The  European  Office  of  Aerospace  Research  and  Development,  London,  England 


DTIC 


1  J.  Blaustcin  Inst,  for  Desert  Research,  Ben-Gurion  University  of  the  Negev,  Sede  Boker  Campus  84990,  Israel  • 

2  Department  of  Civil  Engineering,  Technion-IIT,  Haifa  32000,  Israel. 

3  Department  of  Neurosurgery,  Technion-IIT,  Haifa  32000,  Israel. 

4  Neurosurgeon,  Medical  Center,  POB  133,  Sedgefield  6373,  South  Africa. 

5  School  of  Mechanical  Engineering,  Witwatersrand,  South  Africa. 

6  School  of  Medicine,  Technion-IIT,  Haifa  32000,  Israel.  t-  *  ~  .  '| 

iij  w  '*■  “** 


Unclassified 


SECURII Y  CLASSIFICA1I0N  OF  (HIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704  0188 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4  riRIORMING  OKliANI/AllON  IILPORI  NUMUEH(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Technion  -  Research  and 
Development  Foundation 


lb.  RESTRICTIVE  MARKINGS 


8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 

AFWAL  and  EOARP 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 


RESS  (City, 

AFWL/MLSA  F.0ARD/LRP 

Wright-Patterson  AFB  Box  14 
OH  45433-6533  FPO  New  York  0951 


3  .  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
Distribution  unlimited. 


5.  MONIIOItING  ORGANIZAIION  REI'ORI  NUMIILH(S) 
•  EQftRD  IR  90-09 


7a.  NAME  OF  MONITORING  ORGANIZATION 

European  Office  of  Aerospace  Research  i 
and  Development 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Box  14 

1  FPO  New  York  09510-0200 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFQSR-89-0054 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROJECT 

TASK 

NO. 

NO. 

11.  TITLE  (Include  Security  Classification) 

Cerebral  -  Body  Perfusion  Model 

12.  PERSONAL  AUTHOR(S) 

S.  Sorek,  J.  Bear,  M.  Feinsod 

13a.  TYPE  OF  REPORT 

13b.  TIME  COVERED 

14.  DATE  OF  REPORT  (Year,  Month,  Day) 

15.  PAGE  COUNT 

Annual 

fromNov.88  Tflor.  90 

July  1990 

16.  SUPPLEMENTARY  NOTATION 

» 

COSATI  CODES 


GROUP 


^18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Compartmental  modeling.  Perfusion  flux  and  pressure 
Interaction  between  cerebral,  respiratory  and 
heart  systems;  sensitivity  analysis,  f 


19.  ABSTRACT  ( Continue  on  reverse  if  necessary  arid-ldentify  by  block  number) 

*•  A  lumped  parameter  model  is  developed  to  focus  on  the  dynamic  flow  and  pres¬ 
sure  interactions  between  the  cerebral,  cardiovascular  and  the  respiration 
systems.  The  interrelated  pressures  and  fluxes  are  excited  by  left  cardiac 
pressure,  by  expiration/inspiratiorr  fluxes  and  by  pressure  exerted  on  the 
abdomen.  Sensitivity  analysis  examines  changes  in  pressure  and  flux  at  the 
cerebral  carotid  arteries  and  capillaries  resulting  from  a  sudden  rise  to  an 
upright  position,  changes  in  inhal  e/exhaje  patterns,  presurizing  the  abdomeij, 
changes  in  gravity  acceleration,  changes  In  blood  viscosity  and  heart  fre¬ 
quency.  Animal  experiments  are  performed  to,  validate  model  predictions  and  to 
enable  parameter  estimations,  ) 


20.  DISTRIBUTION/ AVAILA8II IT Y  OF  ABSTRACT 
□  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT, 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr.  Jeff  Higle,  Capt. 


DD  Form  1473,  JUN  86 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

C3  otic  users  Unclassified 


22b.  TELEPHONE  (Include  Area  Code)  22c.  OFFICE  SYMBOL 


Previous  editions  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

Unclassified 


Julius  Silver  Institute  of  Biomedical  Engineering 

Technion,  Israel  Institute  of  Technology,  Haifa,  Israel  32000 


GRANT  AFQSR-89-0054 


CEREBRAL-BODY  PERFUSION  MODEL 


S.  Sorek1,  J.  Bear2,  and  M.  Feinsod3 

in  Collaboration  with 

K.  Allen4,  L.  Bunt5  and  S.  Ben-Haim6 


July  1990 

SCIENTIFIC  REPORT  NO.  4 


Approved  for  public  release,  distribution  unlimited. 


Act  ms  Ion  for 


IIIS  QRA4I 
OIIC  TAB 
Unannounced 
Justification. 


Distribution/ 


Availability  Codas 


jlvail  and/or 
Spatial 


Prepared  for 

The  United  States  Air  Force,  Air  Force  Office  of  Scientific  Research 

and 

The  European  Office  of  Aerospace  Research  and  Development,  London,  England 


1  J.  Blaustcin  Inst,  for  Desert  Research,  Bcn-Gurion  University  of  the  Negev,  Scdc  Bokcr  Campus  84990,  Israel. 

2  Department  of  Civil  Engineering,  Tcchnion-IIT,  Haifa  32000,  Israel. 

3  Department  of  Neurosurgery,  Tcchnion-IIT,  Haifa  32000,  Israel. 

4  Neurosurgeon,  Medical  Center,  POB  133,  Scdgeficld  6573,  South  Africa. 

5  School  of  Mechanical  Engineering,  Witwatcrsrand,  South  Africa. 

6  School  of  Medicine,  Technion-IIT,  Haifa  32000,  Israel. 


Table  of  Contents 


Abstract 

Non-Steady  Compartmental  Model  of  Interactive 
Perfusion  between  Cerebral  and  Body  Systems:  I. Theory 

Abstract 

Introduction 

Model  Topology 

Flow  Equations 

Parameter  Estimation 

Pressure  Waves 

Conclusion 

References 

Non-Steady  Compartmental  Model  of  Interactive  Perfusion  between 
Cerebral  and  Body  Systems:  II.  Sensitivity  Analysis 

Abstract 

Introduction 

Governing  equations 

Sensitivity  analysis 

Conclusion 

References 

The  Effect  of  Positive  Pressure  Ventilation  on  Cordio-Pulmonary  and 
Brain  Dynamics 

Methods 

General  preparation 
Data  acquisition 


Results 


ABSTRACT 


A  lumped  parameter  model  is  developed  to  focus  on  the  dynamic  flow  and  pressure 
interactions  between  the  cerebral,  cardiovascular  and  the  respiration  systems.  The 
interrelated  pressures  and  fluxes  are  excited  by  left  cardiac  pressure,  by 
expiration/inspiration  fluxes  and  by  pressure  exerted  on  the  abdomen.  Sensitivity  analysis 
examines  changes  in  pressure  and  flux  at  the  cerebral  carotid  arteries  and  capillaries 
resulting  from  a  sudden  rise  to  an  upright  position,  changes  in  inhale/exhale  patterns, 
pressurizing  the  abdomen,  changes  in  gravity  acceleration,  changes  in  blood  viscosity  and 
heart  frequency.  Animal  experiments  are  performed  to  validate  model  predictions  and  to 
enable  parameter  estimations. 


NONSTEADY  COMPARTMENTAL  MODEL  OF  INTERACTIVE  PERFUSION 
BETWEEN  CEREBRAL  AND  BODY  SYSTEMS:  I.  THEORY 

S.  Sorek,1  K.  Allen,2  M.  Feinsod,3  S.  Ben  Haim4 
J.  Bear*  and  L.  Bunt6 

ABSTRACT 

A  lumped  parameter  model  is  developed  to  simulate  flow  and  pressure  interaction  between 
the  cerebral  and  the  body  systems.  Its  objective  is  to  study  the  dynamic  interaction  between  the 
cranio-spinal,  body  respiratory  and  the  heart  systems  that  influence  the  brain,  as  well  as  the 
influence  of  conditions  in  the  environment  on  the  body.  By  providing  forecasts  of  departures  from 
a  normal  behavior,  the  model  will  serve  the  following  medical  purposes:  (a)  facilitate  the  under¬ 
standing  of  the  physiology  and  the  mechanisms  that  preserve  the  delicate  brain  in  the  face  of  liv¬ 
ing  stress,  and  (b)  provide  information  for  management  in  deviant  cases. 

The  compartmental  model  consists  of  six  compartments  that  describe  the  cerebral  system, 
and  eight  compartments  that  are  assigned  to  the  body  system.  The  model  also  accounts  for  the 
surrounding  environment  affecting  the  abdominal,  and  respiratory  fluxes.  Altitudes  assigned  to 
the  compartments  introduce  the  effect  of  hydrostatic  pressure. 

Key  Words:  Compartmental  modelling,  interaction  between  cerebral  respiratory  and  heart  sys¬ 
tems:  perfusion  flux  and  pressure:  conductances;  compliances:  environmental 
pressure:  hydrostatic  pressure:  inhale/exhale  fluxes. 

Abbreviated  Title:  Cerebral  and  body  perfusion  model 
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1.  INTRODUCTION 

In  earlier  papers  (Kami  et  al.  1987;  Sorek  et  al.  1988, 1989),  a  compartmental  model  was 
built  to  simulate  the  perfusion  in  the  cerebrovascular  system.  The  mode!  provided  information  on 
pressures  and  fluxes  in  response  to  excitation  in  the  form  of  temporal  flux  and  pressure  changes 
at  the  internal  carotid  artery.  The  influx  to  and  efflux  from  the  cerebrovascular  system  were  intro¬ 
duced,  respectively,  as  external  conditions  imposed  on  the  internal  carotid  artery,  and  the  jugular 
bulb. 

,  In  the  present  work,  this  model  is  extended  to  include  both  the  cerebrovascular  and  the  body 
system  in  a  single  model.  The  objective  is  to  enable  the  study  of  tho  dynamic  interaction  between 
these  two  systems.  In  this  model,  relevant  parts  of  the  brain  and  of  tho  human  body  are 
represented  as  compartments  that  interact  with  each  other,  e.g.,  in  the  form  of  pressure  transmis¬ 
sion  and  exchange  of  fluid.  Each  such  compartment  is  represented  by  lumped,  or  averaged,  pro¬ 
perties  and  state  variables  of  that  organ.  Examples  of  such  lumped  parameters  are  pressure  in  a 
compartment,  inllux/efflux  through  its  boundaries  and  fluid  source/sink  terms. 

The  compartmental  model  is  comprised  of  six  compartments  that  represent  the  cerebral  sys¬ 
tem  and  eight  compartments  assigned  to  the  body  system.  The  body  portion  of  the  compartmen- 
tal  design  involves  only  those  parts  of  the  body  that  are  relevant  to  the  cerebral  system,  i.e., 
interact  grossly  with  the  latter.  Accordingly,  these  parts  include  mainly  the  heart  and  the  respira¬ 
tory  systems. 

The  model  also  accounts  for  the  influence  of  the  environment  surrounding  the  body,  e.g., 
atmospheric  pressure,  pressure  at  high  and  low  altitudes,  high  underwater  pressure,  vacuum  and 
excessive  gravity  acceleration,  such  as  encountered  in  flight  maneuvers.  It  also  takes  into  con¬ 
sideration  differences  in  elevation  between  compartments.  This  means  that  it  takes  into  account 
hydrostatic  pressure. 

Excitations  ol  the  model  are  introduced  through  changes  in  inhale/exhale  rates,activalion 
pressure  of  the  left  cardiac  ventricle  and  outside  (environmental)  pressure  exerted  on  the  abdomi¬ 
nal.  Extreme  intervention,  such  as  resuscitation  and/or  clogging  in  the  heart  system,  can  be  intro¬ 
duced  via  changes  in  the  appropriate  compartmental  conductances  and  compliances. 
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The  diagrammatic  layout  of  the  compartments  is  such  as  to  assist  a  medical  clinician  in  the 
interpretation  of  the  model  image  of  the  human  body,  as  when  examining  an  X-ray  plate.  The 
model  will  guide  clinicians  by  predicting  the  brain-body  perfusion  responses  to  various  excitation, 
thus  enhancing  the  physiological  understanding,  fault  finding  and  consequent  management  in 
deviant  cases  of  this  complex  system. 

2.  MODEL  TOPOLOGY 

.  The  cerebral  section  of  the  model  is  based  on  the  works  of  Kami  et  at.  (1987)  and  Sorek  et 
al.  (1988, 1989),  and  is  comprised  of  the  following  compartments  (Fig.  1): 

Arterial  Cranium  (Ac )  -  Consists  of  four  supply  vessels  through  the  right  and  left  internal 
carotid  and  vertebral  arteries,  with  arteriolar  branches. 

Capillaries  (C)  -  Represents  Ihe  suit  of  brain  capillaries,  choroid  plexus  and  arteriolar/venous 
capillaries. 

Venous  Cranium  (Vc )  -  This  is  a  'lake'  of  blood  confined  by  thin  walls  immersed  in  the  brain 
mass.  It  contains  a  controlled  mechanism  (not  expressed  in  the  mathematics  of  the  model)  to 
influence  quick  and  slow  fluid  movement  in  Ihe  brain  during  stress.  It  comprises  deep  and 
superficial  systems,  normally  freely  anastomising. 

Venous  Sinuses  (S)  -  These  are  encased  in  semi-rigid  walls  to  prevent  collapse  in  all  but 
extreme  conditions  of  compression. 

The  Cerebro-Spinal  Fluid  is  contained  in  two  compartments: 

Ventricles  (F^)  -  The  lour  ventricles  are  treated  as  a  true  compartment  with  its  own  inflow 
and  outflow.  Although,  normally,  the  resistance  to  the  latter  is  low  (as  the  aqueduct  ol  Silvius 
depends  on  a  pulsative  drive  to  maintain  its  patency),  in  disease,  high  resistance  and  even  occlu¬ 
sion  may  occur  and  isolate  part  of  Ihe  ventricular  system.  In  such  cases,  we  introduce  in  the 
model  the  extra  ventricular  and  the  spinal  fluid  as  an  additional  compartment. 

Extra  Ventricular  CSF  +  Spinal  Fluid  (F'EV)  -  The  cisterns  sulci  and  spinal  lluid  compart¬ 
ment  maintains  a  free  communication  with  the  CSF  extra  ventricular  and,  to  a  lesser  extent,  with 
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the  four  main  ventricles. 

Compartments  representing  the  extra-cranial  venous  (systemic)  circulation  are: 

Extra-thecal  venous  plexus  (VET)  -  A  by-pass  of  the  systemic  flow  included  in  the  bony  spi¬ 
nal  portion  of  the  cranial-spinal  compartment.  Normally,  this  compartment  is  much  less  involved 
in  the  circulation  to  the  cerebral  system.  However,  in  special  cases  it  communicates  directly  with 
the  ventricles  through  the  spinal  and  extra  ventricular  fluid.  Therefore,  we  include  it  to  accommo¬ 
date  abnormal  situations. 

.  Brain  Tissue  (B)  -  This  compartment  is  placed  diagrammatically  as  a  central  one,  emphasiz¬ 
ing  its  vulnerability,  especially  due  to  its  location  between  the  venous  compartments.  All  pres¬ 
sures  acting  on  B  must  be  in  such  balance  as  to  preserve  the  tissue  physically,  yet  allow  and  pro¬ 
mote  optimal  macro  and  micro  circulation  of  fluid  transporting  metabolites. 

The  body  section  of  the  model  is  comprised  of  the  following  compartments: 

Respiratory  System  (R)  -  This  compartment  is  activated  by  the  inhale/exhale  flux.  Its 
volume  deforms  through  interaction  with  the  heart  system,  abdominal,  superior  vena  cava  and 
body's  arteries. 

Abdominal  (BD )  -  This  compartment  is  subject  to  volume  changes  initiated  by  the  environ¬ 
ment  pressure,  inferior  vena  cava  and  the  respiratory  system. 

The  heart  system  circulation  receives  an  input  perfusion  from  the  superior  vena  cava  and  a 
feedback  input  from  the  inferior  vena  cava.  It  serves  as  a  circulatory  pump  to  ensure  forward  flow 
of  blood.  Measured  mean  fluxes  are  imposed  on  this  system  so  that  without  respiration  there  still 
is  a  mechanism  forcing  blood  to  How  in  a  specified  circulation.  The  input  flux  is  injected  to  the 
Right  Ventricle  (Rv)  which  mutually  interacts,  by  volume  deformation,  with  Ihe  respiratory  system, 
and  by  flow  with  the  Left  Ventricle  (Lv).  The  left  ventricle  exchanges  volume  deformations  with 
the  respiratory  systems  and  discharges  flow  to  the  body  arteries. 

Body  Arteries  (Ag )  -  This  compartment  receives  its  inflow  from  the  left  heart  ventricle  and 
introduces  the  major  inflow  to  the  cerebral  arteries.  It  interacts,  via  volume  deformation,  with  the 
respiratory  system,  and  discharges  flow  into  the  superior  and  inferior  vena  cava  compartments. 
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Superior  Vena  Cava  (V$q)  -  This  compartment  receives  the  outflow  from  the  cerebral  sys¬ 
tem.  It  also  receives  inflow  from  the  inferior  vena  cava  and  from  body  arteries  and  discharges  it  to 
the  right  ventricle.  Volume  deformation  interaction  exists  with  the  respiratory  system. 

Inferior  Vena  Cava  (V^  )  -  This  compartment  enables  flow  between  body  arteries  and  supe¬ 
rior  vena  cava.  It  interacts  by  volume  deformations  with  the  abdominal,  thus  implicitly  obtaining 
information  from  the  environment  which  is  then  transmitted  as  back  flow  to  the  right  ventricle. 

Next  we  describe  the  resistances  72^  (sRji )  and  compliances  Qy(=CJt  ),  ascribed  to  the 
various  compartments.  Here,  subscript  ij  in  (  )fy  denotes  the  mutual  boundary  of  compartments 
i  and  j. 

RacAb  Carotid  Body  resistance  that  regulates  and  controls  the  inflow  from  the  body  to  the 
cerebral  arteries.  The  control  is  governed  by  the  flux,  QA ,  to  the  cerebral  system. 
CacB  Compliance  between  the  brain  tissue  and  cranial  arteries  compartments. 

CacFzv  Compliance  that  attenuates  the  arterial  pulse  transmitted  by  large  vessels  traversing 
the  basal,  cisternal,  sulcal  spaces. 

RacC  1116  sum  °* the  cranial  arteries,  capillary  and  choroid  plexus  resistances.  The  resis¬ 
tance  associated  with  the  capillary  is  under  autoregulatory  control.  This  control 
between  the  arteriolar  and  capillary  vessels,  attenuates  the  systole  artery  pulse.  The 
autoregulatory  effect  will  not  be  reflected  in  the  mathematics  of  the  model. 

Rqvc  Resistance  of  the  Arteriolar/Venous  capillary,  accounting  for  the  pressure  drop 
observed  between  them. 

RCB  Resistance  of  the  Blood-Brain  barrier  (between  the  capillary  and  the  brain  tissue). 

RcFiy  Endothelial  resistance  of  the  Blood-CSF  barrier.  It  describes  the  choroid  plexus  and 

ependymal  secretion. 

CcF,y  Compliance  manifesting  the  arterial  pulse  transmitted  to  the  ventricle  CSF  by  the 
choroid  plexus  and  the  extra-cellular  fluid. 
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RVcB 

RVCS 


RF„B 

CBF„ 


Resistance  representing  the  Blood-Brain  barrier  (involved  in  cerebral  oedema). 

A  resistance  comprised  of  two  gates.  First,  the  outlet  from  the  deep  venous  circula¬ 
tion  into  the  straight  sinus  via  the  great  cerebral  vein  of  Galen.  According  to  Le  Gros 
Clark,  an  auto-regulatory  control  may  be  associated  with  this  outflow. 

The  second  component  is  the  superior  cortical  veinous  outflow  into  the  sagittal  sinus. 
An  observable  pressure  drop  exists  across  this  resistance  under  normal  conditions. 

Ependymal  resistance  of  the  CSF-Brain  barrier. 

Multifrequency  pulsed  compliance,  transmitted  rhythmically  by  the  brain  to  the  ventri¬ 
cle  for  axial  drive  of  CSF. 


RF/vFgv  Resistance  that  exists  only  in  the  case  of  strong  stenosis,  between  the  ventricles  and 
the  extra  ventricular  compartments. 

rFhvS  Resistance  manifesting  villous  tufts  secretion  into  the  venous  sinus. 

CsFgy  A  low  Compliance  across  the  semi-rigid  sinus  walls. 

RFf/Vsc  Resistance  that  manifests  the  slow  secretion  area  around  the  spinal  root  sleeves. 

CfkyVe T  Compliance  between  the  Intradural  spinal  fluid  and  the  systemic  venous  flux,  across 

the  dura  and  the  extra-thecal  venous  plexus.  This  mediates  postural,  respiratory, 
abdominal  and  other  body  fluctuations,  while  setting  up  the  intracranial  pressure 
level. 


RVgrVx  Resistance  to  secondary  flow  communication  between  the  extra-thecal  venous 
plexus  and  the  superior  vena  cava. 
rVxab  Resistance  to  inflow  from  body  arteries. 

Cvxr  Compliance  created  by  the  respiratory  venous  flux.  This  compliance  is  controlled  by 
the  involuntary/voluntary  variations  in  the  respiratory  rhythms. 

Ry K rv  Resistance  between  the  superior  vena  cava  and  the  heart  system  via  the  right  ventri¬ 


cle 
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Ryx  yK  Back  flow  resistance  into  the  superior  vena  cava  from  the  inferior  one. 

CrrV’  CrLv 

These  are  the  compliances  between  the  respiratory  and  the  heart  system,  comprised 
of  the  right  ventricle,  and  the  left  ventricle,  respectively. 

Crah  Compliance  of  the  respiratory  and  body  arteries. 

Crbd  Compliance  between  the  respiratory  and  the  abdominal.  It  transmits  the  environmen¬ 
tal  effects  to  the  abdominal. 

RlvAh  Resistance  to  flow  discharged  from  the  left  ventricle  (leaving  the  heart  system)  into 
the  body  arteries. 

Rabvk  Back  How  resistance,  from  the  body  arteries  into  the  inferior  vena  cava,  to  become  a 
flow  feedback  to  the  heart  system. 

C 'bdvk  Compliance  between  the  abdominal  and  the  inferior  vena  cava,  transmitting  the 

environmental  effect  to  the  abdominal. 

Cbda  Compliance  due  to  direct  communication  between  the  surrounding  environment  and 
the  abdominal. 

Next,  we  write  the  perfusion  equations  for  the  entire  cerebral-body  compartmental  setup. 

3.  FLOW  EQUATIONS 

Let  us  consider  a  single  incompressible  fluid  phase  that  approximately  represents  all  the 
relevant  fluids  in  the  brain-body  system. 

In  writing  the  fluid  flow,  or  balance  equations,  wo  consider  mass  and  momentum  balances 
for  each  compartment  of  the  cerebral-body  compartment  system. 

Essentially,  each  balance  equation  states  that  the  temporal  rate  of  increase  of  either  the  fluid 
mass,  or  its  momentum  in  a  compartment,  is  equal  to  the  amount  of  net  influx  of  that  quantity 
through  the  compartment’s  boundaries  plus  the  external  sources  (e.g.  injection  or  ejection  of  fluid) 
within  the  compartment.  For  a  constant  density  fluid,  the  mass  balance  reduces  to  a  volume 
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balance. 


The  fluid's  volume  balance  equation  in  compartment  n  takes  the  form 


in 


where  Vn(t)  denotes  the  volume  of  fluid  in  compartment  n  at  time  t  (equal  to  the 
compartment's  volume),  Qn  is  the  source  associated  with  compartment  n ,  and  gm-  denotes  the 
flux  flowing  out  of  compartment  n ,  through  the  boundary  of  the  compartment,  to  its  adjacent  i 
compartment.  The  non  rigidity  of  the  compartment's  boundaries  (*walls)  is  expressed  by  a  com¬ 
pliance  factor,  Cnj ,  defined  by 


(2) 


where pnj(=pn-pj )  denotes  the  pressure  difference  between  compartments  n  and  j,  on  both 
sides  of  their  common  wall. 

For  low  Reynolds  number  flow,  the  momentum  balance  of  a  fluid  moving  through  a  capillary 
tube  can  be  shown  to  reduce  to  an  equation  that  expresses  linear  proportionality  between  flux  and 
driving  force  (Sorek  et  at.  1988).  The  latter  is  composed  of  a  pressure  gradient  and  a  gravity 
term.  Here  we  assume  that  a  similar  expression  governs  the  flux  between  adjacent  compart¬ 
ments.  Hence 

Qni  =  Zni  (Pni  +  Y#m  )  =  W  I  „/  (3) 


where  (zh )  |  nJ-  (s  zni  ( hn  -hi )),  Zni  is  the  conductance  associated  with  the  flow  wall  bolwoen 
compartments  n  and  i ,  //m-  (=//„  )  is  the  altitude  difference  between  compartments  n  and 

Pni 

i,  y  is  the  fluid's  specific  weight,  /tm-  ( =  +  Hn{ )  denotes  the  piezometric  head  difference 

between  compartments  n  and  i  and  zm-  ( =  7 )  is  called  the  hydraulic  conductivity  factor 
between  compartments  n  andi. 
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Upon  substituting  (2)  and  (3)  into  (1),  we  obtain 

im\ni+nCp)\nj=Qn  (4) 

»'  j 

where  Pnj(=~ (Pn~Pj ) =  Pn~Pj )  denotes  the  temporal  rate  o(  increase  in  the  ditference  in 

pressure  between  compartments  n  and  its  adjacent  one,  j. 

Under  certain  abnormal  conditions,  a  number  of  compartments  become  active.  These  are 
the  extra  ventricular  and  spinal  fluid  (F^y)  and  extra-thecal  venous  plexus  (Vet).  In  the 
mathematical  model,  this  fact  is  expressed  by  introducing  a  parameter  Xmn  which  can  be  set 
either  to  X^  =  0  for  normal  conditions,  or  to  Xnm  =  1  for  abnormal  ones. 

During  the  pressure  cycle,  a  situation  may  oocur  in  which  the  flow  can  change  ils  direction 
through  two  adjacent  compartments  as  a  result  of  change  in  their  pressure  difference.  To  prevent 
this,  we  introduce  a  value  term  denoted  by  T]nm  which,  accordingly,  may  have  the  values 
T)nm  =  0  for  normal  flow  direction,  or  Tinm  =  1  for  abnormal  one. 

In  previous  reports  (Sorek  et  al.  1988,  1989)  the  flow  problem  of  the  cerebral  system  was 
solved  using  arterial  pulsatile  flow  as  a  boundary  condition.  The  pulsatile  nature  of  the  flow 
(arterial  and  all  other  vessels  for  that  matter)  is  now  transmitted  via  the  cyclic  behavior  of  the 
piezometric  head. 

In  view  of  (4),  we  write  the  following  compartmental  fluid  balance  equations: 

Ac:  0  =  —{zh ) |4„4C  +  (zh ) “*■  {Pp ) \a<:B  "f" (XCp )  (5.1) 

C:  0  =  -(z/i)j4cc  +{z^)!co  +{zh)\cvc  +  {zh)\cF„  +  [Cp)\cFiY  (5-2) 

Vc:  0  =  -{zA)|cvc-(2A)|Bvc +  (z/»)lvcs +(Cp)|veu  (5.3) 


B  :  0  as  -{zh  )lcc-(2*)l  f„d  +{zh)\  bvc+{Cp  )  I  uvc+(Cp )  I  bf„+(Cp  ) !  BAc 


(5.4) 
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S:  0  =  (z/t) | sv*;  -(2Mlvcs  —  ) I FzyS  +  {Cp)\sfiv  (5.5) 

Fjv'  0  =  -{zh  )  I  CF^zfl  )  i  FyB  +(kzh  )  I  FvFtv+iCP  )  I  F„B  MPp  ) !  /VvC  (5.6) 

Fgv".  0  =  {zh)\F„S  +(^)lFfyV*r-(^)lF^riv  +  (^P)li!,«vS  + 

(5.7) 

+  (kCp )  I  FrvAc  +  {XCp )  I  f„vff 

Vet ’  0  =  (teA )  |  v^Vac  +  (^Qp )  I  v,rF*v  (5.8) 

Vsc :  0  =  -{z/l )  |  SVsc  ~  (zfl )  I  V/cV*e  “  izfl )  I  Ai  V«  +  (^  )  I  V»c/?y  “ 

(5.9) 

+  (teA )  I  VrTVx  -  (teh )  I  FYv Vac  +  (Qj )  |  vmR 

VIC:  0  =  -{zh ) |/,vK  +  {zh)\ vKvx  +  (Cp ) ! VkBo  (5.10) 

Ab:  Q  =  {zh)\A,Ac-{zh)\LvA'l  +(z/i)Utv«: +(^)li4,V/c+(Qp)U»«  (5.11) 

Rv:  0  =  -{2/i)|v*;fiy  +  (^)lfiviv  +  (<^P)!flvH  (5.12) 

Lv:  o  =  -{n?A)|fiv£,y  +  (n*/i)l£^,  +  (C!p)Iz,v«  (5.13) 

R  '•  Qr  =  {Cp )  |  RVK+{Cp )  I  RRv+(Cp )  I  Rr+{Cp )  I  RLy  +  (Qp )  I  RDo+(PP )  I RA,  (5.14) 

Bd\  0  =  \Cp)  I  DnR  +{Cp)\DDvlc+{Cp)\BnA  •  (5.15) 


For  the  compartment  representing  the  environment,  we  can  write  the  relation 
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(Cp  )\ab„+Qr~^  (6.1) 

If  the  environment  is  of  infinite  extent  (i.e.  atmosphere),  we  may  write 

Pa=  0.  (6.2) 

We  assume  a  Monro-Kellie  postulate  according  to  which  the  (almost)  rigidity  of  the  skull  dic¬ 
tates  that 

.  ^Qsdt  =jQAdt  (7) 

where  T  is  the  time  period. 

A  stenosis  in  the  passage  between  the  Fjy  and  Fgy  compartments,  initiates  a  build-up  of 
the  compliance  between  the  Fjy  and  Aq  compartments.  This,  in  return,  indicates  that 

I-F/vFkv  ~  FevAc ^  •  (8.1) 

Under  normal  conditions,  with  free  communication  between  Fjy  and  FFy  compartments 
0),  both  the  main  ventricle  and  the  extra  ventricular  will  merge  into  one  compart¬ 
ment  (the  F  compartment),  resulting  in 

hFfV  =  hFjy  m  hF  with  1  =  0  .  (8.2) 

Note  that  in  (9)  we  have  assumed  HF/V  =HFtysHF,  because  only  minor  differences  exist 
between  their  values. 

To  prevent  a  build  up  of  pressure  in  compartment  Ly  higher  than  that  of  Ry,  we  introduce 
a  value  term  denoted  by  Similarly  denotes  the  value  between  compartments  .A# 

andLy.  The  values  are,  therefore, 
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0  :  hLv<  hAv 

l\RvLv  -  (9-1) 

1  '<  hiy  2  hRv 
0  I  hLy  ^  hAa 

T\LyAB  ~  (9-2) 

1  ;  hLv  <  hAg 


(11.1) 
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0 

0 

0 

0 


0 

0 


nCg,Qn  +  CvmgCB,Ci,ghi,  +  CyKg  Cfcjf  Cfl„^  hA  ] 


_1_ 

c 


Qn  +  Cv*B,  Cg,  If  Cu ghif  +  Cy^jj,  Cg  C/j,^  hA  j 


(11.2) 


^UAt^AtU^lr  +  £•  [CajR^&Q/?  +  ^AtHCBtCuR^U  +  C^,nCg,[iCi]B^h/i] 

^n.u^nrLr^u  +  ~j=  [Cn,HcBmQR  +  CrtR  Cb»Ci*r  hi*  +  Crtr  Cb,r  ^bba  J 


The  non  zero  elements  of  the  z  matrix  are: 


10,10  =  CrvR  [1  — 
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The  sum  of  the  abdominal  compliances,  CBj>,  and  that  of  the  respitory,  CR  ,  are  given  by 


CB0  =  CBdR  +  CBdVk  +  CB0A 

(11.5) 

cr  =  CRVsc  +  Cmv  +  CjtLy  +  CRBd  +  CRab 

(11.6) 

Also, 

C  =  CBd  Cr  -  C§dR  .  (1 1 .7) 

Equation  (10)  is  to  be  solved  for  the  unknown  h  vector,  subject  to  an  initial  condition  of  the 

form 

h.(t=o) = ho  •  (12) 

Upon  solving  (10),  together  with  (12),  we  obtain  the  solution  for  the  temporal  rates  of  change 

♦  • 

in  piezometric  head,  hRo  and  hR ,  in  the  forms 
h Bo  =  -=  [CBD  r  Qr  +Cbd  RCRAB  ^Ab  +Cr  Cbd  v,c+C«  Cbda  + 

(13) 

+CbdR  CRVxhvSc+CBDR  Crrv ^Ry+Cso R  CRLV^LV  ] 


hR  ~  i  IC0O QR+CBn CRA„hAB+CBDR  CBd Vl(hvK+ 


CBd CRVschv9C^BD CRRvhRv+CBDRCBDA +  CBDCRLvhLv  )  . 


(14) 


The  compartmental  balance  equation,  written  in  the  global  compact  form  (10),  Involves  con¬ 
ductivities  and  compliances  expressed  by  the  matrices  z  and  C,  respectively.  In  order  to  solve 
for  the  piezometric  head  vector,  h ,  we  need  to  know  the  z  and  C  values.  We  shall,  next,  discuss 
an  inverse  method  for  estimating  these  parameters. 
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4.  PARAMETER  ESTIMATION 

To  predict  the  pressure  and  flux  response  of  the  model  to  external  changes,  the  values  of 
the  parameters  C  and  z  must  be  known.  In  order  to  estimate  them,  we  need  measured  values  of 
pressure  in  the  various  compartments  at  a  sufficient  number  of  points  in  time.  Typical  phasic 
pressures  are  depicted  in  Figs.  2-11. 

In  the  present  work,  we  assume  that  all  conductivities  and  compliances  have  the  form  of  a 
time  step  function  i.e.,  they  change  from  one  constant  value  to  another  constant  value,  due  to 
abnormal  situations  such  as  disease.  Because  C  and  z  are  constant  during  long  time  intervals, 
and  4  (0  >s  assumed  cyclic,  by  taking  a  temporal  average  of  (10),  i.e.,  integration  over  a  period  of 
time  divided  by  the  period,  we  obtain 

(15) 

where  h*  and  Q*  denote  mean  values  of  the  piezometric  head  and  of  the  source  term,  respec¬ 
tively.  Since  h*  and  Q*  are  known  averaged  values,  in  order  to  solve  for  z  values,  we  rewrite 
(15)  in  the  form 

£V=£‘  (16) 

where  f)  denotes  the  head  difference  between  two  communicating  compartments. 

In  the  cerebral  portion  of  the  model,  we  face  ten  conductivities  against  six  balance  equa¬ 
tions.  However,  by  virtue  of  the  Monro-Kellie  doctrine  which  assumes  absolute  rigidity  of  the  cra¬ 
nial  vault,  actually  only  five  equations  of  the  six  are  independent.  Thus  for  this  redundancy,  five 
additional  conditions  are  needed.  One  of  them  is  the  mean  influx  to  the  CSF  (F)  compartment 
Qh  o  .3ml/min.  (Cutalar,  1968).  Two  more  conditions  are  stipulated  by  physiological  data. 
These  conditions  determine  the  scalar  coefficients  (Sorek  et  al.  1987)  a  and  P  defined  as 

a=fM_  (a  >  0)  (17) 

ZFD 


and 
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P  =  _i±_  (0  ^  p  <  oo)  (18) 

ZFB 

where  a  («  10000)  indicates  the  ratio  of  the  vein-venous  sinous  to  the  cerebrospinal  fluid-brain 
barrier,  and  P(  =  1 0-3)  is  the  ratio  of  the  biood-brain  barrier  to  the  cerebrospinal  fluid-brain  bar¬ 
rier. 

The  mean  arterial  influx,  QA ,  which  according  to  (3)  equals  the  mean  venous  sinous  efflux, 
is  given  by 

'  Qa=Qs  =  Q*  =750ml./min. 

Because  Qp,  Qa  and  Q$  are  a-priori  known,  and  the  mean  head  values 
hAcc<  hcF<  ^AflAc  and  hsv3C  are  a,so  known> the  following  conductivities  are  immediately 
obtained. 

Q*  Q* 

zAcAb  =  7? -  zAcC  =  TT”  (19.1) 

** AgAc  ■ AcC 


(19.2) 


Also,  for  the  given  configuration  of  the  heart  system  (Fig.  1),  we  impose  the  direction  and 
magnitude  of  the  same  mean  flux  through  conductances  of  the  cardiac  system  ensuring  a  circula¬ 
tory  pump  effect.  Hence,  such  a  condition,  e.g.,  Q[  =  5900ml./min  through  the  boundary 
between  compartments Ly  andAB ,  will  yield 

QL 

zLvAg  =  -  (20-1) 

hLvAa 

QL 

zVscRv  ,  * 

hVxRv 


(20.2) 
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We  assume  that  the  mean  influx  to  the  inferior  vena  cava  compartment  is  a  fraction,  8 
(0  <  8  £  1 ),  of  the  efflux  QL  from  the  left  heart  ventricle.  Therefore,  we  obtain 

SQL 

zvlcvsc  =  7 ; -  (20.3) 

hvKvx 


SQL 

2a„vk=t: —  (20.4) 

hABVlc 


(1-5  )QL-Q* 

zAbVsc  = "  *  (20.5) 

hABV9C 


Hence,  in  view  of  (17)  to  (20),  equation  (16)  is  constructed  of  the  following  matrix  and  vec¬ 


tors 


r 


h Br  ~  P^Cfl 

0 

0 

hBVc 

P  hCB 

hCVc 

0 

0 

K FB 

0 

h-FS 

0 

o^vcs 

~hCVc 

0 

-hLv 

(21.1) 


ZV  =  lzFB  <  zCVc>  ZFS  >  ZBVC? 


(21.2) 


Q'  ={0,Q*-QL"Qf>  Of  (21.3) 

Upon  solving  the  inverse  problem  for  the  conductivities,  we  move  to  the  estimation  of  the 
compliances. 

Given  information  of  the  simultaneous  pulse  wave  recordings,  p^j,  andp^j,  in  the  different 
compartments  and  at  various  times  t  *  ,k  =1 ,2 equation  (10)  (or  (5)  and  (6))  now  yields  a 
set  of  k  relations  for  the  compliances,  for  times  t  * 
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(Zh)k  lA'Ac-izh^c 

[2hk)\AcC-(2hk)\cB-(2hk)  \cvc-{zhk)\CF 

(zhk)\cvc  +  {zhk)\BVc-zhk  I  vcs 
{zhk)  I  vcs  +  (zhk) \FS  -  (zh k )  | $ya 

-(zhk)\vc$+(zhk)\FB+(zhk)  \A ERc-(zhk)\AcCHzhk)  | cB«?hh)\cvc 

-lzhh)\sv9C  Hzhk) \AbAc 

(zhk)\AaAc  +  (2^*)li4,Vjc  +  (zhk)\v,cV9c  ~  (2^*,l VxRv 
{zhk)\  Vb  vK~(zhk)\  VKVX 
(zhk)\vscRv-(i]zhk)\RvLv 

(nzhk)\RvLv-(nzhk)\LvA8 

~l?hk)\ABvK-(zhk)  \AaVac  - {zhk)\AgAc  +  {r]zhk)\LvAn 
Q*+(zhk)  I vKVsc 
Qr  +zv,oABhABvK 


Note  that  the  z  values  in  (23.3)  are  already  known. 

By  the  Gauss-Markov  theorem,  the  C  values  are  derived  as  an  assortment  of  the  set  of 
informations  through  all  K  time  observations  (Sorek  et  al.  1987b),  namely 


C  =  (kt  n)-1  kTB 


where, 

P 

h2 


n  = 


K*  13 


(24) 


(25.1) 


(23.3) 


B=[R\R2 . RK],lK 


Since  the  hk  submatrices  of  it  are  main  diagonalized,  we  can  evaluate  the  C  values 


directly. 


CAcB  = 


2 

k  77  “ 

£  Kb? 

*=  1 


Crv  — 


£  hcFm  * )  UcC  -  (*A* )  I CB  -  W * )  I CV0  -  * )  I  Cf] 

*=1 

_  — — 

£  K? 

k= 1 


£  Kb K*A*)lcvc  +  (*A*)layc - (*A*)|Vcs] 

CVcB  *“77  - 

2  (avcb)2 


2  AsF[^*)lvcS+(^*)lFS“^*)|SVa] 
r  _  *=1 _ 

CSF~  K  -i  9 

2  feO2 

*=1 


Cnp  = 


2^  -  4b  i(zh  k  )  I  VCS  ~{zh  * )  I FB  ~(zh  k )  I  AnAc+(zh * )  I AcC~(zh  * )  I  CB ~(zfl * )  I  CVc] 

x  7TZ  '  _ 

2  (4b)2 

*= i 


(26.5) 
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(26.11) 
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K  ..  . 

X  hRBD  l Qr  +  ZVScV,MkV9c 

CBdR  =  — - F— -  (26.12) 

X  rbd  )2 

*=1 

X  ^Bda[Qr  +{zhk)  \abvk 

CBdA  =  i=! - ^ -  (26.13) 

X  (*b0a)2 

A=1 

Thus,  we  have  completed  writing  the  numerical  algorithms  for  estimating  the  z  and  C 
values.  The  mean  and  temporal  values  of  the  fluxes,  pressure  and  pressure  rates  are  derived 
from  phasic  changes  as  depicted  in  Figs.  2-1 1 . 

Note  that  for  the  body  system,  the  comparlmentat  pressure  and  temporal  rate  ol  pressure 
change  ride  on  the  appropriate  respiratory  wave.  However,  for  the  cerebral  system,  this  effect  is 
attenuated. 

We  assume  a  respiratory  pressure  wave  of  the  form 

Pa  =  10  sin  Wflf  (27) 

2jt 

where  0)^  =  is  the  inspiration/respiration  frequency,  and  TR  ( =  3.4sec.)  is  taken  as  the 

1R 

respiratory  cycle  time. 

We  use  the  following  altitude  values: 


HBd  *  110  cm. 

Hv„  «  140  cm. 
y  sc 

Hyc  « 163  cm. 

Hr  » 135  cm. 

//y^a  126  cm. 

Hb  « 163  cm. 

HRv  ■  130  cm. 

HAb  *  116cm. 

Hp  » 163  cm. 

HP  *  130  cm. 

HAc  m  160  cm. 

Hs  » 163  cm. 

Hlv  *  130  cm. 

Hq  ■  1 63  cm. 
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corresponding  to  a  standing  position.  Also,  we  choose  the  following  values  for  pressure  and  rates 
of  pressure  changes: 

pc  =  30 mm.Hg.  ;  pc  =  0 
Ps  =  2 mm.Hg.  ;  p,  =  0 
Pv,c  =  5mm. Hg.  ;  pyK  =  0  . 

This  concludes  the  calibration  scheme  for  the  brain-body  compartmental  model,  with  con¬ 
ductivities  and  compliances  that  remain  constant  during  long  time  intervals. 

5.  PRESSURE  WAVES 

Once  the  linear  model  described  by  (10}  has  been  calibrated,  it  is  possible  to  obtain  solu¬ 
tions  for  the  pressure  waves  in  all  brain-body  compartments.  In  incremental  form,  the  solution  of 
the  differential  matrix  equation  (10)  for  h  is 

=  expMK^o-J'1^)]  +J"1£(0  •  (28) 

Here,  A/  ( =  tm+^-tm)  denotes  a  time  difference  between  the  frontier  time  level,  tm+*, 
and  the  backtime.  t m .  The  matrix  t  reads 

x  =  A*  C_1|  (29) 

The  rational  expansion  of  exp(-x ),  leads  to  the  following  formula 
/-(1-6)x 

exp(-t)  =  (Osesi)  (30) 

in  which  /  is  the  unit  matrix.  Substituting  (29)  and  (30)  into  (28),  yields 

4(1*4.)  =  U+6'  r'LHi-e)t  Kio-r'S#)!  +I",210  •  <31> 

The  coeflicient  0  controls  the  type  of  numerical  solution  scheme  evolving  in  time.  When 
0  =  0,  we  have  an  explicit  scheme  0  =  1,  an  implicit  scheme,  and  0  <  0  <  1  is  the  mixed 
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scheme  case. 

Thus,  with  the  choice  of  0,  the  head  waves  in  the  various  compartments  are  calculated 
according  to  (10).  In  view  of  (3)  and  (31),  we  can  also  evaluate  the  compartmental  pressures 

£{<+a<)  =  )]  +  2_1fit(t)-S} Y  (32) 

This  completes  the  modelling  of  the  brain-body  perfusion  pressure,  as  excited  by  environ¬ 
mental  pressure  and  due  to  inspiration/expiration  flux. 


6.  CONCLUSION 

A  compartmental  model  for  perfusion  pressure  interplay  between  the  cerebral,  respiratory 
and  heart  systems  is  presented.  The  perfusion  pressure  response  is  initiated  by  compartments 
altitude,  environmental  pressure  and  inspiration/expiration  flux. 

Flow  between  compartments  is  governed  by  compartmental  step  function  related  conductivi¬ 
ties  and  compliances  that  take  the  forms  of  time  step  functions. 

Once  the  calibrated  model  is  available,  it  will  be  used  for  the  prediction  of  phasic  pressures 
associated  with  perfusion,  and  the  influence  of  various  situation  such  as: 

I:  Flight  maneuvers,  subjecting  the  body  to  acceleration  forces  of  several  times  the  gravity 
acceleration  [g ),  in  the  direction  of  head  to  foot. 

The  idea  here  is  to  increase  the  gravity  head  (which  is  part  of  the  hydrodynamic  head)  gradually, 
in  order  to  check  its  effect  on  the  resulting  perfusion  fluxes,  especially  the  influx  to  the  arterial 
cranium  and  those  in  the  cerebral  system.  Parameters  influencing  the  deterioration  of  cerebral 
flux  curves,  manifesting  the  so-called  ‘blackout’,  can  also  be  investigated.  For  example,  the  effect 
can  be  studied  of  different  inhale/exhale  and/or  different  relative  positioning  between  head-body, 
in  overcoming  the  characteristics  of  phasic  pressure  curves  resulting  from  the  excessive  g 
values. 

II:  High/low  altitudes  and  vacuum,  generating  extreme  environmental  cases. 

Such  situations  affect  the  pattern  of  inhale/exhale  (amplitude  and  periodicity).  These  affects  on 
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the  phasic  pressure,  can  be  studied,  taking  into  account  possible  changes  in  fluid  density.  In  addi¬ 
tion,  at  high  altitudes,  the  fluid’s  viscosity  varies  and  this,  in  turn,  affects  the  conductance 
(inversely). 

Ill:  Different  inhale/exhale  maneuvers. 

We  intend  to  examine  the  response  of  phasic  pressures  in  the  various  model  compartments,  for 
example,  to  artificial  inhale/exhale,  as  applied  in  resuscitation. 

IV:  Different  relative  positioning  of  the  Brain-Body  as  in  physical  exercises.  We  intend  to 
ch$ck  the  influx,  for  example,  into  the  common  carotid  artery,  due  to  tilting  the  body  relative  to  the 
head.  Another  example  is  the.case  of  standing  upside  down. 

V:  Different  degrees  of  occluding  the  heart  and  body  conductances. 

This  relates,  for  example,  to  resuscitation  procedures  in  which  the  abdominal  and  the  respiratory 
systems  are  brought  under  cyclic  pressure  variations.  In  our  model,  this  may  cause  changes  in 
the  conductances  and  compliances  of  the  body  model. 

The  developed  model  can  be  used  to  guide  the  neurosurgeon  in  the  interpretation  of  possi¬ 
ble  consequences  of  management  methods  applied  to  abnormal  cases.  Following  are  some 
examples  of  management  maneuvers  that  can  be  interpreted  by  numerical  experiments  of  the 
model. 

(a)  Correction  of  a  fault  if  recognizable  and  accessible,  e.g.,  tumours,  blockage  of  flow  chan¬ 
nels.  This  is  the  ideal  case. 

(b)  Offsetting  of  destructive  tensions  that,  usually,  would  be  of  short  duration,  e.g.,  reducing 
excess  pressure  usually  done  by  shunting  the  CSF  volume. 

(c)  Correcting  mechanisms  that  fail  to  confine  or  maintain  the  system  in  optimal  yielding  against 
repeated  stress  in  the  long  term,  e.g.,  introducing  a  gas  bubble  in  a  sac  into  the  large  hydro¬ 
cephalic  ventricle  in  order  to  hold  down  the  pressure  peaks  and  create  better  Pressuro- 
Volume-Time  relationship. 

In  practice,  the  management  of  (b)  is  on  the  CSF  bulk,  where  as  for  (c)  it  also  involves  the 
blood  volume  as  an  indirect  means  of  promoting  optimum  flux. 
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The  model  can  provide  information  for  fault-finding,  target  identification  and  monitoring  the 
results  of  rational  management. 
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Figure  1:  Cerebral  -  Body  compartmental  design. 
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Figure  2:  During  one  cycle  of  breathing 
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TYPICAL  PRESSURE  IN  B 


Figure  3:  During  one  cycle  of  breathing 
(after  Hamit  et  al.) 
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Figure  4:  During  one  cycle  of  breathing 
(after  Hamit  et  al.) 
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Figure  5:  During  one  cycle  of  breathing 
(after  Hamit  et  a l.) 
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TYPICAL  PRESSURE  IN  Ab 


Figure  6:  During  one  cycle  of  breathing 
(after  Shepherd  et  at) 
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Figure  7:  During  one  cycle  of  breathing 
(after  Wiggers) 
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Figure  8:  During  one  cycle  of  breathing 
(after  Vernon  et  al) 
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Figure  1 0: During  one  cycle  of  breathing 
( after  Vernon  et  al.) 
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NONSTEADY  COMPARTMENTAL  MODEL  OF  INTERACTIVE  PERFUSION 
BETWEEN  CEREBRAL  AND  BODY  SYSTEMS  II.  SENSITIVITY  ANALYSIS 

S.  Sorekl,  O.  Lichtenstein^,  J.  Bear^,  M.  Feinsod4  and  S.  Ben-Haim^ 

Abstract 

A  sensitivity  analysis  is  performed  in  a  lumped  parameter  model,  focusing  on  the 
dynamic  flow  and  pressure  interactions  between  the  cerebral,  cardiovascular  and  the  pulmonary 
systems. 

The  interrelated  pressures  and  fluxes  are  excited  by  left  cardiac  pressure,  by 
expiratory/inspiratory  fluxes  and  by  pressure  exerted  on  the  abdominal.  The  analysis  examines 
changes  in  pressure  and  flux  at  the  cerebral  carotid  arteries  and  capillaries  resulting  from  a 
sudden  rise  of  the  body  to  an  upright  position,  changes  in  inhale/exhale  flux  patterns, 
pressurizing  the  abdomen,  changes  in  gravity  acceleration,  changes  in  blood  viscosity  and 
changes  in  heart  frequency.  Results  are  consistent  with  available  clinical  information  and  the 
outcome  may  be  explained,  qualitatively,  on  the  basis  of  understanding  the  physiological 
mechanisms  and  on  clinical  observations. 

Key  Words:  Compartmental  modelling;  perfusion  flux  and  pressure;  cerebral,  cardiac  and 
respiratory  systems;  conductances;  compliances;  gravity  acceleration;  cardiac  and  abdomen 
pressure  amplitude  and  frequency;  carotid  and  capillary  to  brain  tissue  influx. 
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A.  Introduction 

In  the  first  paper  (Sorek  et  al.  1990)  a  model  was  presented,  in  which  the  brain-body 
system  was  conceptualized  as  a  compartmental  configuration,  representing  lumped  vessels  and 
ventricles  of  similar  physiological  functions  in  transmitting  perfusion  flow.  Physical  laws, 
expressing  fluid's  mass  and  momentum  balances,  were  then  formulated  and  combined  to  form 
a  set  of  equations  that  relate  pressures  to  fluxes,  using  time  step  functions  of  conductances  and 
compliances.  The  same  concepts,  when  applied  to  modelling  brain  mechanics  (Sorek  et  al., 
1988a)  suggested  a  possibility  of  excessive  fluid  accumulation  in  the  ventricles,  due  to 
occluding  perfusion  between  capillaries  and  brain  tissue  (i.e.  blood-brain  barrier).  Actual  brain 
surgery  of  an  elderly  patient  suffering  from  Normal  Pressure  Hydrochephalus  (NPH),  verified 
this  prediction  (Sorek  et  al.,  1988b).  Conductances  and  compliances  were  then  estimated  by  an 
inverse  method  (Sorek  et  al,  1989b).  With  these  estimated  values,  temporal  variations  of 
pressure  in  the  various  compartments  were  derived  by  solving  the  model,  subject  to  a  given 
pressure  cycle  at  the  cerebral  artery  (Sorek  et  al.,  1988c).  The  simulated  temporal  pressure 
variations  of  the  Jugular  Bulb  and  that  of  the  Venous  Sinus,  demonstrated  suction  intervals 
during  the  pressure  cycle.  It  was  demonstrated  that  increasing  drainage  through  a  shunt 
implanted  in  the  ventricles  may  cause  intensified  deterioration  of  the  characteristic  pressure 
curve  of  the  ventricular  CSF.  These  features  are  consistent  with  clinical  observations  and, 
together  with  previous  results,  justified  the  use  of  the  compartmental  brain  mechanics  model, 
CBMM  (Sorek  et  al.,  1988a). 

Next,  a  multi  compartmental  brain-body  (MCBBM)  (Sorek  et  al.,  1990)  was 
constructed  to  account  for  various  situations,  such  as: 

A)  Flight  maneuvers  subjecting  the  body  to  dynamic  forces  of  several  g  (gravity 
acceleration)  magnitudes  in  the  direction  of  head  to  foot. 

B)  High/low  altitudes  and  suction,  generating  extreme  environmental  cases. 

C)  Different  respiratory  maneuvers. 

D)  Different  relative  positioning  between  the  brain  and  the  body,  as  in  physical  exercises. 

E)  Different  patterns  of  pressuring  the  abdomen. 
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F)  Simultaneous  application  of  a  combination  of  the  above  maneuvers,  as  in  the  case  of 
resuscitation. 

Measurements  of  the  actual  changes  in  hemodynamics  during  these  maneuvers  is 
reported  in  the  literature.  Blood  pressure  changes  during  orthostatic  testing  (Wolthuis  et  al., 
1979;  Fulco  et  al.,  1985);  alteration  in  heart  rate  following  linear  accelerations  (Lauritzen  et  al., 
1987);  perfusion  pressures  as  well  as  intracranial  pressure  distribution  at  different  head 
elevations  (Rosner  and  Coley,  1980);  changes  in  the  venous  pressures  during  orthostatic 
testing  (Jonson  and  Rundcrantz,  1969)  as  well  as  alteration  in  cardiac  output,  total  peripheral 
resistance  following  head-down  tilt  (London  et  al.,  1983)  were  all  reported  and  are  a  short 
listing  of  recent  experimental  literature.  Most  of  these  reports  cannot  delineate  the  direct  effect 
of  the  maneuver  studied  on  the  cardio  vascular  system  from  the  compensatory,  reflex  response 
brought  about  by  the  antonomic  nervous  system.  Our  approach  may  supply  an  analysis 
framework  with  which  the  direct  effect  of  the  studied  perfurbation  can  be  assessed. 

The  compartmental  setup,  as  described  in  the  first  paper,  is  comprised  of  six  cerebral 
compartments  and  eight  compartments  assigned  to  the  body  system  (Fig.  1).  Excitation  of  the 
model  is  imposed  via  respiratory  fluxes,  by  activation  pressure  of  the  left  cardiac  ventricle  and 
by  external  pressure  exerted  on  the  abdomen. 

Let  us  first  briefly  discuss  the  theoretical  base  for  constructing  the  flow  equations  for 
the  brain-body  compartmental  assemblage.  The  sensitivity  analysis  will  then  follow,  with  the 
objective  of  assessing  the  influence  of  the  aforementioned  maneuvers,  subject  to  the  prescribed 
excitations  driving  the  model. 

B.  Governing  equations 

Practically,  within  the  range  of  regular  brain-body  functioning,  we  may  regard  its 
interconnecting  fluid  phases  (i.e.  blood  and  Cerbro  Spinal  Fluid,  CSF)  as  being 
incompressible,  not  a  function  of  constituents'  concentrations  and  independent  of  temperature. 

Therefore,  for  each  compartment,  we  consider  the  mass  and  momentum  balance 
expressions  for  a  single  fictitious  fluid  phase  of  constant  density. 
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Each  such  balance  equation  states  that  the  temporal  increase  of  an  extensive  quantity  of 
the  fluid,  here  mass  or  momentum,  within  a  compartment  is  equal  to  the  net  influx  through  the 
compartment's  boundaries,  plus  the  resultant  sum  of  the  external  sources  of  the  considered 
quantity. 

For  a  constant  density  fluid,  and  for  a  compartment  with  a  finite  time  dependent 
volume,  the  fluid  mass  balance  reduces  to  a  balance  equation  in  the  form 

dVn 

"dt~ + X  qni  *  X  =  * Wn  M 

j  i 

where  Vn(t)  denotes  fluid's  volume  within  compartment  n  at  time  t,  qnj  denotes  the  efflux  from 
compartment  n  to  the  adjacent  j  one,  through  their  common  boundary,  and  qin  denotes  influx  to 
compartment  n  from  compartment  i,  through  their  common  boundary;  Qn  denotes  a  source  due 
to  the  injection  fluid  into  compartment  n  and  wn  denotes  a  sink  ejection  fluid  out  of 
compartment  n.  A  storage  argument  is  introduced  as  a  result  of  the  non-rigidity  of  the 
compartment  boundaries.  This  is  expressed  by  a  time  step  compliance  factor,  Cnl,  say  at  the 
mutual  boundary  between  compartment  n  and  1,  defined  by 

r  i  dVn  ns 

Cnl "  d(Pnl)  (2) 

where  Pnl(=  Pn  -  Pi)  denotes  the  pressure  difference  between  compartments  n  and  1,  on  both 
sides  of  their  common  boundary. 

For  a  viscous  fluid  moving  slowly  through  a  capillary  tube,  the  compartmental 
momentum  balance  reduces  to  a  linear  proportionality  between  the  flux  crossing  a 
compartmental  boundary  and  the  pressure  and  elevation  difference  on  both  sides  of  the 
boundary.  Hence,  for  the  flux  leaving  compartment  n  into  compartment  j  through  their 
common  boundary  nj,  we  write 


Qnj  =  Znj  (Png  +  pgHnj),  (3.1) 
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where  Znj  denotes  a  conductance  factor  between  compartment  n  and  j;  p  denotes  the  fluid 
density;  g  denotes  the  gravity  acceleration,  and  Hnj  (=  Hn  -  Hj)  denotes  the  elevation 
difference  between  compartments  n  and  j.  Similarly,  for  the  flux  entering  compartment  n  from 
its  adjacent  one,  i,  through  their  common  boundary,  in,  we  write. 

<lin  =  Zin  (Pin  +  pgHin).  (3.2) 

The  lumped  conductance  Znj  (or  Z\n)  is  assumed  constant  for  long  time  periods. 

Actually,  expressions  (3)  relate  to  the  differential  form  of  the  fluid's  linear  momentum 
balance  equation  known  as  Darcy's  law,  which  is  used  in  a  continuum  approach  to  flow 
through  a  porous  medium.  These  equations  may  be  interpreted  as  the  integration  of  Darcy's 
law  over  a  volume  surrounding  a  common  boundary  of  two  compartments  interacting  by  flow 
(Adar  and  Sorek,  1989a).  Hence,  like  the  permeability  appearing  in  Darcy’s  law  so  is  the 
conductance  a  function  inversely  proportional  to  fluid’s  viscosity  (p). 

Upon  substituting  (2)  and  (3)  into  (1),  we  obtain  a  flow  equation  for  each  compartment 
n,  in  the  form 


M 

o 

HI* 

+ 

— i 

(Xznj  +  XZinl 

Pn 

fX'C  1  A 

■  LCnl  dt H 

h  X  ZnjPj  +  XZinPi 

i  i 

■U  1  J 

_ 

u 

j  i  ) 

=  Qn  -  Wn  +  pg  f£ZinHin  -  ]£ZnjHnjN 

.  »  j 


.  (4) 


Note  that  a  change  in  g  levels  introduces  a  modification  of  the  driving  force  (i.e.  the 
l.h.s.  of  (4)).  A  change  in  viscosity  will  produce  an  inverse  change  in  all  of  the  conductance 
terms  on  both  sides  of  (4), 

Given  an  initial  setup  of  pressures  in  the  various  compartments  we  solve  the  set  of 
equations  (4)  for  the  temporal  pressure  values,  subject  to  prescribed  time  dependent  values  of 


respiratory  flux,  external  pressure  exerted  on  the  abdomen  and  activation  pressure  on  the  left 
ventricle. 


C.  Sensitivity  analysis 

The  objective  of  this  sensitivity  analysis  is  to  assess  the  clinical  outcome  when 
subjecting  the  brain-body  system  to  changes  in  various  control  variables  (e.g.  changing  the 
relative  positioning  between  body  and  head),  or  changing  the  amplitude  and  frequency  of  the 
driving  excitations  (e.g.  cardiac  activation  pressure,  or  external  pressure  applied  on  the 
abdomen).  We  focus  on  the  resulting  pressure  at  the  carotid  arteries  and  the  cerebral  capillaries, 
as  well  as  on  the  influx  into  the  arteries  and  from  the  capillaries  into  the  brain  tissue.  In  doing 
so,  we  may  understand  and  predict  clinical  syndromes  based  on  brain-body  flow 
communications,  such  as  the  "black  out"  phenomenon  experience  by  pilots  in  excessive  g 
flights,  or  various  resuscitation  maneuvers. 

Actually,  in  order  to  solve  for  the  pressures  in  (4),  from  which  we  can  then  obtain  the 
fluxes  by  employing  (3),  we  need  to  know  the  true  values  of  the  conductances  and  the 
compliances.  To  evaluate  these  values,  we  should  resort  to  some  inverse  process,  such  as  the 
one  described  by  Sorek  et  al.  (1989b). 

The  estimation  of  these  parameters  will  be  reported  in  the  future.  Instead,  we  assume 
here,  based  on  average  fluxes  and  pressures,  a  set  of  values  for  the  parameters.  These  values 
are  so  selected  as  to  meet  the  prescribed  flow  direction  between  the  compartments.  Since  the 
compliances  and  conductances  are  assumed  to  be  constant  within  long  time  periods,  we  may 
check  sensitivity  expressed  in  terms  of  deviation  from  normal  (average)  variations  of  the 
selected  state  variable.  Denoting  a  as  the  time  dependent  state  variable  (flux,  or  pressure),  we 

express  the  normal  value  of  a  («*Yby 


a  -  «ave 
«ave 


(5) 


where  aave  denotes  the  average  value  of  a  over  a  time  period,  Tr,  defined  by 
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Tr 


«ave  = 


(6) 


We  start  our  sensitivity  analysis  by  asking  whether  we  can  verify  the  model  against  results 
known  qualitatively. 

Obviously  the  model  does  not  include  control  mechanisms.  For  example,  in  a  real 
physiological  system  we  may  sometimes  observe  an  immediate  influence  due  to  neoronical, 
biochemical  and  mechanical  regulation  on,  say,  the  pressure  difference.  The  model,  however, 
describes  flow  interrelations  and,  as  such,  its  predictions  are  important  in  showing  the  result 
that  may  take  place  when  control  is  malfunctioning. 

The  sensitivity  to  a  sudden  rise  of  the  body  to  an  upright  position  is  shown  in  Figs.  2 
and  3.  Note  the  general  trend  of  reduction  in  pressure  in  the  body  arteries,  Ab  (Fig.  2a), 
cerebral  arteries  Ac  (Fig.  2b),  cerebral  capillaries,  C  (Fig.  3a)  and  in  the  cerebral  ventricles,  F 
(Fig.  3b),  similar  to  results  measured  by  Rosner  et  al.  (1980).  This  serves  as  a  reasonable 
validation  of  the  reliability  of  the  developed  model. 

The  sensitivity  to  changes  in  respiratory  flux  is  shown  in  Fig.  4.  Note  that  a  50% 
change  in  amplitude  yields  only  slight  deviations  at  the  end  of  the  diastole  in  the  influx  from 
body  to  brain  (Fig.  4a),  as  well  as  in  the  flux  from  capillaries  to  brain  tissue  (Fig.  4b). 
Hyperventilation  results  in  the  reduction  of  intracranial  pressure,  expressed  in  reducing  the 
volume  associated  with  changes  of  high  inspiratory  volume  and  frequency,  but  not  in 
substantial  flux  changes.  Therefore,  Fig.  4  describes  a  valid  situation  within  the  accepted 
metabolism. 

The  sensitivity  of  flow  in  the  carotid  arteries  and  that  between  capillaries  to  brain  tissue, 
to  changes  in  cardiac  driving  pressure  is  shown  in  Figs.  5a  and  5b,  respectively.  Note  that  a 
higher  cardiac  pressure  may  cause  an  amplification  of  brain  tissue  impedance,  which,  in  return, 
will  yield  a  lower  capillary  to  brain  flux  (Fig.  5b). 

The  sensitivity  to  external  pressure  applied  to  the  abdomen,  is  described  in  Fig.  6.  The 
importance  of  synchronization  with  the  driving  heart  cycle  is  demonstrated  by  comparing  the 
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normal  flow  at  the  carotid  arteries  (Fig.  6a)  and  the  one  from  capillaries  to  brain  tissue  (Fig. 
6b).  Note  that  although  a  higher  external  pressure  is  exerted,  the  sudden  relief  of  the  abdomen, 
creates  a  suction  at  some  portions  of  the  cycle  and  a  flow  buildup  at  others.  This  suggests  an 
intensive  pressure  rise  followed  by  a  slow  relief,  accommodated  by  synchronization  with  heart 
cycle. 

The  sensitivity  to  changes  in  gravity  acceleration  and  its  impact  on  pressure  and  flow,  is 
shown  in  Figs.  7,  8  and  9,  respectively.  Two  extreme  maneuvers  are  exemplified.  The  first 
simulates  an  uplift  with  a  4g  acceleration.  The  second  represents  the  simulation  of  going 
downward  at  an  acceleration  of  0.5g.  In  the  first  case,  the  body  will  experience  a  resultant 
acceleration  of  5g  directed  from  foot  to  head,  while  in  the  second  one,  the  resultant  acceleration 
of  0.5g  will  prevail.  Note  that  in  the  5g  case  (Figs.  7  and  8),  a  pressure  buildup  to  and  over  the 
normal  peak  during  the  systole  is  indicated.  However,  a  steep  fall  is  exhibited  during  the 
diastole.  This  indicates  a  suction  action  which  may  reduce  the  oxygen  level  at  that  instant  The 
o.5g  maneuver  lowers  the  pressures  accordingly,  but  causes  a  moaerate  fall  from  the  systole 
peak  to  diastole.  A  similar  behavior  is  demonstrated  with  regard  to  flow  into  the  carotid  arteries 
and  to  the  one  between  capillaries  and  brain  tissue  (Fig.  9).  Note  the  suction  towards  the  end  of 
the  diastole  period  (Fig.  9a),  associated  with  a  steeper  (compared  to  the  0.5g  curve)  fall  in 
flow.  Fig.  9b,  showing  the  5g  case,  strongly  suggests  a  possible,  so-called,  "luxury 
perfusion",  in  which  natural  anatomic  shunts  (collaterals)  are  opened  to  by-pass  sudden 
amplification  of  capillaries  to  brain  tissue  flow.  In  such  a  case,  the  tissue  does  not  receive  a 
sufficient  supply  of  nutrients  and  the  outcome  may  be  the  so-called  "blackout"  phenomenon. 

The  sensitivity  to  change  in  blood  viscosity,  affecting  carotid  and  capillaries  to  brain 
tissue  flow,  is  described  in  Figs.  10a  and  10b,  respectively.  Lowering  blood  viscosity,  which 
may  take  place  at  high  altitudes,  will  reduce  the  flow  into  the  carotid  arteries  (Fig.  10a),  but 
will  increase  the  capillaries  to  brain  tissue  flow  (Fig.  10b).  This,  in  return,  may  deplete  the 
intake  of  oxygen  by  the  tissue,  which  may  result  in  a  feeling  of  dizziness  experienced  at  high 
altitudes. 

Actually,  Fig.  10a  represents  a  prediction  that  may  contradict  clinical  observations.  We 
may  explain  this  as  an  artifact  of  not  using  correct  conductances.  Lowering  viscosity  increases 
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the  conductances  and  may  thus  result  in  lowering  the  influx  into  the  carotid  arteries  at  the 
expense  of  influx  through  other  boundaries  into  surrounding  compartments. 

The  sensitivity  of  the  carotid  influx  to  simultaneous  changes  in  the  frequency  of  cardiac 
activation  pressure  and  magnitudes  of  the  gravity  acceleration,  is  shown  in  Fig.  11.  Note  that 
the  relative  flow  surfaces  depicted  in  Figs.  11  and  16  are  the  flow  mean  value  (i.e.,  time 
integral  over  respiration  cycle  of  temporal  flow  divided  by  respiration  time  interval)  for  each 
combination  of  g  level  and  heart  frequency  and  abdominal  pulse  frequence,  respectively.  Two 
cross-sections  of  the  relative  carotid  flow  surface  of  Fig.  11  are  described  in  Figs.  12.  It  is 
demonstrated  (Fig.  12a)  that  when  experiencing  an  excess  of  5.5g,  lowering  of  heart  frequency 
will  result  in  gaining  higher  flow  into  the  carotid  arteries.  On  the  other  hand,  when  in  values 
lower  than  the  normal  g  (g%  <  0),  this  flow  remains  constant  up  to  approximately  50%  of  the 
regular  heart  frequency.  Two  sections,  at  high  and  low  heart  frequencies  (fig.  12b),  predict  a 
turning  point  (at  about  lg  value)  of  a  higher  flow  for  a  higher  g  magnitudes  and  comparably 
lower  flow  rate  for  g  values  smaller  than  the  normal  one.  This  is  more  pronounced  for  lower 
frequency  values. 

The  sensitivity  of  the  carotid  influx  to  simultaneous  changes  in  the  amplitude  of  the 
pulse  pressure  exerted  on  the  abdomen  and  the  g  magnitudes,  is  described  in  Fig.  13.  Figure 
14a  demonstrates  that  for  excessive  g  values,  the  carotid  flow  will  increase,  while  increasing 
the  external  pressure  amplitude.  This  phenomenon  is  reversed  when  considering  g  values 
lower  than  normal.  Fig.  14b  exhibits  the  same  phenomenon  as  that  described  in  Fig.  12b. 
However,  note  the  decrease  in  carotid  flow  beyond  twice  the  regular  g  value.  The  carotid  flow 
surface  as  a  function  of  the  frequency  of  the  outside  pressure  applied  to  the  abdomen  and 
change  of  g  levels,  is  depicted  in  Fig.  15.  Figure  16a  describes  (in  the  case  of  excessive  g 
values)  a  trend  of  increasing  the  carotid  flow  when  the  outside  pressure  frequency  is  lowered. 
Note  also  the  support  in  raising  the  flow  when  this  frequency  meets  the  cardiac  frequency. 
Figure  16b  is  similar  to  Fig.  14b. 
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D.  Conclusion 

A  sensitivity  analysis  was  carried  out  on  a  compartmental  unsteady  perfusion  model 

that  describes  the  brain-body  flow  system. 

The  sensitivity  analysis  was  applied  to  various  maneuvers  in  order  to  check  the 

influence  of  changing  control  variables  and  excitations  on  pressure  and  flow  levels  at  the 

carotid  arteries  and  transmission  between  capillaries  and  brain  tissue. 

The  results  indicate  that: 

-  A  sudden  rise  causes  the  lowering  of  pressure  at  the  cerebral  arteries,  capillaries  and 
ventricles. 

-  Changes  in  inhale/exhale  pattern  is  consistent  with  the  Hyperventilation  phenomenon  that 
causes  the  reduction  of  inhale  volume  and  does  not  change  significantly  vascular  influxes. 

-  Changes  in  cardiac  activation  pressure  produce  a  reduction  of  flow  between  capillaries  and 
brain  tissue. 

-  When  subjecting  the  abdomen  to  outside  pressure,  it  is  important  to  synchronize  the  cycle 
with  that  of  the  heart 

-  Excessive  magnitudes  of  g  cause  an  increase  of  influx  from  capillaries  to  brain  tissue.  This 
may  be  interpreted  as  the  luxury  perfusion  situation  which  bypasses  the  flow  with  the 
nutrients  it  carries  (through  newly  opened  collaterals)  and  result  in  a  "blackout" 
phenomenon. 

-  In  low  viscosity  situation  (e.g.  when  at  very  high  altitude)  the  flow  between  capillaries  to 
brain  tissue  increases.  This  (although  not  to  the  same  degree  as  in  excessive  g  magnitudes) 
may  cause  a  sense  of  dizziness  due  to  the  depletion  of  some  oxygen,  following  the  above 
mentioned  mechanism. 

-  Simultaneous  changes  of  heart  frequency  and  gravity  acceleration  predict  that  at  excessive  g 
levels  the  lowering  of  heart  frequency  will  yield  the  rise  of  influx  into  the  carotid  arteries. 

-  Simultaneous  change  of  amplitude  of  pressurizing  the  abdomen  and  changes  of  g 
magnitudes  yield  the  rise  of  carotid  influx  when  increasing  pressure  amplitude,  at  excessive 
g  levels. 
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-  Simultaneous  change  of  the  frequency  of  pressurizing  the  abdomen  and  changes  of  g 
magnitudes  predicts  the  increase  of  the  carotid  flow  upon  lower  pressure  frequency. 

Now  that  we  have  assessed  the  behavior  of  the  brain-body  model  under  various 
conditions,  the  next  step  will  be  to  evaluate  values  for  the  various  conductances  and 
compliances.  This  will  be  reported,  separately,  in  the  future. 
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RELATIVE  CAROTID  FLOW 
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Figure:  13 
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THE  EFFECT  OF  POSITIVE  PRESSURE  VENTILATION  ON 
CARDIO  PULMONARY  AND  BRAIN  DYNAMICS 


In  order  to  obtain  experimental  validation  to  the  model  predictions  and  to  enable  parameter 
estimation  we  have  performed  animal  experiments.  The  objective  of  the  experimental  setup  was 
to  enable  us  to  record  simultaneous  pressure  and  flow  signals  from  different  sites  in  the 
cardiovascular,  pulmonary  and  neural  system.  Recordings  were  performed  under  different 
pressure  and  volume  perturbation  introduced  into  the  thorax.  We  studied  12  large  dogs. 

Methods 

General  Preparation 

Mongrel  dogs  (18-25kg)  were  anesthetized  with  pentobarbital  (60mg/kg  IV),  using 
supplemental  doses  (30mg)  as  dictated  by  presence  of  corneal  reflex  or  pressor  response  to 
surgical  incision.  Animals  were  intubated  and  ventilated  (Harvard  large  animal  respirator)  to 
maintain  physiological  PO2  (80-110)  and  PCO2  (35-45).  Bolus  injections  of  sodium 
bicarbonate  (2mEq/ml)  will  be  used  to  maintain  arterial  pH  between  7.37  and  7.43.  Rectal 
temperature  was  monitored  and  maintained  between  37  and  40°  C,  using  a  heating  pad  as 
needed. 

Cannulae  were  placed  in  the  right  femoral  artery  and  vein  for  measuring  arterial  pressure  in  the 
descending  aorta  and  administration  of  fluids  and  pharmacological  agents,  respectively. 
Another  set  of  cannulae  were  placed  in  the  left  femoral  artery  and  vein  for  measuring  arterial 
blood  gases  and  inferior  vena  cava  pressures,  respectively  (Fig.  1). 

Via  a  ventral  incision  in  the  neck  the  common  carotid  artery  and  the  internal  jugular  vein  were 
isolated.  Electromagnetic  flow  probes  of  appropriate  size  were  placed  on  these  vessels  for 
measuring  blood  flow. 

Through  an  incision  above  the  parietal  bone  a  burr  hole  (2  mm)  was  drilled,  keeping  the  dura 
matter  intact.  After  hemostassis  was  performed  using  bone  wax,  a  Millar  micro  tip  needle 
pressure  transducer  (3F)  was  introduced  into  the  brain  tissue  for  measuring  tissue  pressure 
(Fig.  2).  Through  a  puncture  of  the  cisterna  magna  another  miliar  micro  mannometer  was 
introduced  for  measuring  cisternal  pressure. 
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Air  way  pressure  was  measured  using  Microswitch  pressure  transducer  connected  to  a  PE200 
tube  placed  near  the  intratracheal  tube  opening.  Airflow  through  the  airways  was  measured 
using  a  Gould  pneumotachograph  connected  at  the  intratracheal  tube  connection  to  the  positive 
pressure  ventilator  (Fig.  3). 

Pressure-Volume  perturbation  were  introduced  by  a  computer  controlled  positive  pressure 
ventilator.  This  ventilator  was  built  using  a  bellow  connected  to  a  stepper  motor.  The  maximal 
displacement  of  the  bellow  as  well  as  the  acceleration  and  deceleration  of  its  movements  were 
externally  controlled  by  a  computer.  An  algorithm  enabled  on-line  modification  of  two  of  these 
parameters  (maximal  displacement  and  acceleration).  The  deceleration  of  the  bellow  as  well  as 
the  respiratory  rate  were  calculated  by  the  controlling  computer  to  ensure  constant  minute 
ventilation.  For  evaluation  the  actual  minute  ventilation  was  measured  continuously  by 
integration  of  the  pneumotachograph  inspiratory  signal. 

Data  acquisition 

Analog  signals  were  fed  into  a  12  bit  A/D  converter,  and  sampled  at  1  Khz.  Data  were  stored  in 
a  digital  form  for  off  line  processing. 

Results 

12  dogs  were  studied.  A  typical  tracing  of  the  pressure  and  flow  signals  during  normal  positive 
pressure  ventilation  is  depicted  in  Fig.  4. 

It  is  evident  that  the  pressure-volume  perturbation  introduced  into  the  thorax  had  a  significant 
and  observable  effect  on  each  of  the  pressure  and  flow  signals  recorded. 

In  general  the  pressure  responses  were  similar  to  a  low  pass  filter  operation  to  the  applied 
airway  pressure,  with  a  varying  degree  of  delay.  Arterial  blood  flow  into  the  brain  was 
pulsatile  but  was  superimposed  on  the  lower  frequency  of  the  respiratory  perturbation.  The 
venous  flow  was  continuous  but  varied  periodically  following  the  respiratory  perturbation. 
Figure  5  depicts  five  typical  recordings  of  the  pressure  flow  time  curves  during  ventilation  with 
(a)  trapezoid,  (b)  saw  tooth,  (c)  triangle,  (d)  low  inspiratory  to  expiratory  ratio  as  well  as  (e) 
high  inspiratory  to  expiratory  ratio. 
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Figure  Legends: 

Fig.  1  Experimental  setup  (panel  a)  together  with  a  schematic  representation  of  sensors 

position  in  the  corresponding  theoretical  model  map  (panel  b).  1  -  Aortic  pressure,  2  - 
Carotid  flow,  3  -  Intraventricular  pressure,  4  -  Intracerebral  tissue  pressure,  5  -  Jugular 
vein  flow,  6  -  Superior  vena  cava  pressure,  7  -  Intrathoracic  pressure  (pleural),  8  - 
Inferior  vena  cava  pressure,  9  -  Intrabdominal  pressure. 

Fig.  2  Parietal  bone  hole  prepared  for  inserting  intracranial  pressure  transducer. 

Fig.  3  Pressure  and  flow  transducers  connected  to  the  airway  opening. 

Fig.  4  A  typical  tracing  of  the  pressure  and  flow  signals  during  normal  positive  pressure 
ventilation. 

Fig.  5  Typical  recording  of  the  pressure  flow  time  curves  during  ventilation  with  (a)  trapezoid, 

(b)  saw  tooth,  (c)  triangle,  (d)  low  inspiratory  to  expiratory  ratio  as  well  as  (e)  high 
inspiratory  to  expiratory  ratio. 
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Figure  1:  Cerebral  -  Body  compartmental  design. 


FIGURE  4  NORMAL  (FLOW)  VENTILATION 
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FIGURE  5b  SAW  TOOTH  (FLOW)  VENTILATION 
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FIGURE  5d  LOW I:E  RATIO 
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